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Neutrino fluxes from proto-neutron stars with and without quarks are studied. Observable differ- 
ences become apparent after 10-20 s of evolution. Sufficiently massive stars containing negatively- 
charged, strongly interacting, particles collapse to black holes during the first minute of evolution. 
Since the neutrino flux vanishes when a black hole forms, this is the most obvious signal that quarks 
(or other types of strange matter) have appeared. The metastability timescales for stars with quarks 
are intermediate between those containing hyperons and kaon condensates. 
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A proto-neutron star (PNS) is born following the grav- 
itational collapse of the core of a massive star, in conjunc- 
tion with a successful supernova explosion. During the 
first tens of seconds of evolution, nearly all (~ 99%) of the 
remnant's binding energy is radiated away in neutrinos 
of all flavor [Q. The ^—luminosities and the evolutionary 
timescale are controlled by several factors, such as the to- 
tal mass of the PNS and the v— opacity at supranuclear 
density, which depends on the composition and equa- 
tion of state (EOS). One of the chief objectives in mod- 
eling PNSs is to infer their internal compositions from 
v— signals detected from future supernovae like SuperK, 
SNO and others under consideration, including UNO || . 

In their landmark paper, Collins and Perry J3j noted 
that the superdense matter in neutron star cores might 
consist of weakly interacting quarks rather than of 
hadrons, due to the asymptotic freedom of QCD. The ap- 
pearance of quarks causes a softening of the EOS which 
leads to a reduction of the maximum mass and radius 
||. In addition, quarks would alter v— emissivities and 
thereby influence the surface temperature of a neutron 
star H during the hundreds of thousands or millions of 
years that they might remain observable with such in- 
struments as HST, Chandra and XMM. Quarks would 
also alter the spin-down rates of neutron stars |j] . 

Many calculations of dense matter predict the appear- 
ance of other kinds of exotic matter in addition to quarks: 
for example, hyperons or a Bose (pion, kaon) condensate 
(cf. |Q and references therein) . An important question is 
whether or not v observations from a supernova could re- 
veal the presence of such exotic matter, and further could 
unambiguously point to the appearance of quarks. The 
detection of quarks in neutron stars would go a long way 
toward the delineation of QCD at finite baryon density 
which would be complementary to current Relativistic 
Heavy Ion Collider experiments, which largely address 
the finite temperature, but baryon-poor regime. 

An important consequence of the existence of exotic 
matter in neutron stars (in whatever form, as long as 
it contains a negatively charged component), is that a 
sufficiently massive PNS becomes metastable ||. After 
a delay of up to 100 s, depending upon which component 
appears, a metastable PNS collapses into a black hole 



]9|, |10yil| . The collapse to a black hole proceeds on a free- 
fall timescale of less than a millisecond JT2J , much shorter 
than v— diffusion times, and the neutrinos still trapped 
in the inner regions cannot escape. Such an event should 
be straightforward to observe as an abrupt cessation of 
v— flux when the instability is triggered M]. 

The evolution of the PNS in the so-called Kelvin- 
Helmholtz phase, during which the remnant changes from 
a hot, v— trapped, and lepton-rich object to a cold and 
v— free star, occurs in near-hydrostatic equilibrium. The 
v— matter interaction timescales are much smaller than 
the dynamical timescale of PNS evolution, which is of the 
order of seconds. Thus, until neutrinos enter the semi- 
transparent region at the edge of the star, they remain 
close to thermal equilibrium with matter, and may be 
treated in the diffusion approximation. 

In this Letter we provide a benchmark calculation with 
quarks by solving the general relativistic ^—transport 
and hydrostatic equations (as in JlOj] ) , and then com- 
pare our results with those of our previous work jlO 11 
in which other compositions were studied. In addition, 
we assess the prospects of observing PNS metastability 
and its subsequent collapse to a black hole through ob- 
servations in current and planned detectors. 

The essential microphysical ingredients in our study 
are the EOS of dense matter and its associated 
v — opacity. We begin by considering two generic compo- 
sitions: charge-neutral, beta equilibrated matter contain- 
ing (i) nucleons only (np) and (ii) nucleons with quark 
matter (npQ). In the npQ case, a mixed phase of baryons 
and quarks (pure quark matter exists only for very large 
baryon densities, except for extreme choices of model pa- 
rameters) is constructed by satisfying Gibbs' phase rules 
for mechanical, chemical and thermal equilibrium | |l4) . 
The EOS of baryonic matter is calculated using a field- 
theoretic model at the mean field level |15|. The results 
reported with this EOS are quite general, as we verified 
by alternatively using a potential model approach 
The quark matter EOS is calculated using a MIT bag- 
like model (similar results are obtained with the Nambu- 
Jona-Lasinio quark model). The details of the EOS may 



be found in Ref. Il6|. We use ^—opacities |17 
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tent with the EOS. When quarks appear, the v absorp- 
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tion and scattering cross sections dramatically decrease, 
the precise reduction being sensitive to the thermody- 
namic conditions in the mixed phase Jig]. 



mass. The appearance of quarks is accompanied by an 
increase in Y v because of the depletion of electrons; for 
the largest masses, the increase is very large. 
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FIG. 1. Evolutions of the central baryon density tib, v con- 
centration Y Vy quark volume fraction \ an d temperature T 
for different baryonic masses Mb- Solid lines correspond to 
stable stars; stars with larger masses are metastable (dashed 
lines). Diamonds indicate when quarks appear at the star's 
center, and asterisks denote when metastable stars become 
gravitationally unstable. 

Fig. |l| shows the evolutions of some thermodynamic 
quantites at the center of npQ stars of various, fixed, 
baryonic masses (Mb)- hi the absence of accretion, Mb 
remains constant during the evolution, while the gravita- 
tional mass Mq decreases. With the EOS utilized, stars 
with Mb J$ 1-1 Mq do not contain quarks and those with 
Mb ^ 1.7 M are metastable. The latter value is ~ 0.05 
Mq larger than the maximum mass for cold, catalyzed 
npQ matter, because the maximum mass of hot v— free 
npQ matter is this much less than that of hot i/— trapped 
matter |fT8fl . Generally, due to the high lepton (y) con- 
tent initially present in the PNS, the electron chemical 
potential at the center is too large for quarks to exist. For 
sufficiently massive stars, quarks eventually appear after 
a certain amount of v loss occurs. For the Mb = 1.6 
Mq star, for example, quarks appear after about 15 s 
(indicated by a diamond). Thereafter, the star's central 
density increases for a further 15-20 s, until a new sta- 
tionary state with a quark-hadron mixed phase core is 
reached (for stable stars) or for Mb ^ 1-7 M Q , instabil- 
ity occurs (indicated by asterisks). It is interesting that 
for this EOS the lifetimes for all masses are restricted to 
the range 10-30 s, and slowly decrease with increasing 
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FIG. 2. A comparison of i? e count rates expected in SuperK 
from a PNS containing either np or npQ matter. The left 
panel shows times less than 10 s, while the right panel shows 
times greater than 10 s. 

To point out the major differences one might observe 
between the np and npQ cases, we have estimated the v e 
count rate in the SuperK detector in Figure [| For this 
estimate, we assumed the total ^—luminosity from a PNS 
at 8.5 kpc distance, corresponding to a Galactic super- 
nova, was equally divided among the six v species. The 
v — energy spectra were taken to be Fermi-Dirac with zero 
chemical potential and a temperature corresponding to 
the matter temperature in the PNS where the ^—optical 
depth was approximately unity. Only the dominant re- 
action, v e absorption on protons, was included. (For de- 
tails, see 0.) It is difficult to discern much difference in 
the early (t < 10 s) count rates from np and npQ stars. 
For stars with Mb < 1.7 M , this is because quarks have 
not yet appeared. For more massive stars, the fact that 
neutrinos are strongly trapped inhibits any discrimina- 
tory signal from reaching the surface before this time. 
The signals at later times (t > 25 s), however, are sub- 
stantially larger for the npQ case, due to the decrease in 
v— opacity of npQ matter and the increased binding en- 
ergy of npQ stars. Most importantly, the v— signal from 
metastable npQ stars halts abruptly when the instability 
occurs. Qualitatively, these features are also found for 
npH and npK stars ■ 

We compare ^—signals observable with different de- 
tectors in Fig. ^, which displays the v— light curves as 
a function of Mb for npQ stars. The two upper shaded 
bands correspond to estimated SN 1987A (50 kpc dis- 
tance) detection limits with KII and 1MB, and the lower 
bands correspond to estimated detection limits in SNO, 
SuperK, and UNO, for a Galactic supernova (8.5 kpc dis- 
tance). The detection limits have been set to an v e count 
rate dNJdt — 0.2 Hz |ll| with the same assumptions as 
in Fig. g. The general rise with time of the detector lim- 
its is chiefly due to the steady decrease in the v e average 
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energy. It is possible that these limits are too conserva- 
tive and could be lowered with identifiable backgrounds 
and knowledge of the direction of the signal. The width 
of the bands represents the uncertainty in the v e average 
energy due to the flux-limited diffusion approximation. 
We conclude that it should be possible to distinguish be- 
tween stable and metastable stars, since the luminosities 
when metastability is reached are always above conser- 
vative detection limits. 
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FIG. 3. The total ^—luminosity for npQ stars of various 
baryon masses. Shaded bands illustrate the limiting luminosi- 
ties corresponding to count rates of 0.2 Hz for the indicated 
supernovae in some detectors. 

The drop in v— luminosity for stable stars is associ- 
ated with the end of the Kelvin-Hclmholtz epoch when 
the PNS is becoming optically thin. This portion of the 
v— light curve is approximate due to the breakdown of 
the diffusion approximation. It is an apparent coinci- 
dence that this occurs simultaneously with the collapse 
of the lower mass metastable stars. 

Our choice of bag constant, B — 150 MeV fm -3 , in 
conjunction with the baryonic EOS we used, was moti- 
vated to maximize the extent of the quark matter phase 
in a cold neutron star, and was limited by the necessity 
of producing a maximum mass cold star in line with ac- 
curate observational constraints (Mq = 1.444 M [[l9||). 
Increasing B, or employing an alternative quark EOS 
that otherwise produces a larger maximum mass, delays 
the appearance of quarks and raises the metastability 
window to larger stellar masses. Necessarily, this results 
in an increased timescale for metastability for a given 



mass, and a lower z^— luminosity when metastability oc- 
curs. Fig. H shows the relation between time to instabil- 
ity and Mb for the original case (B — 150 MeV fm -3 , 
thick solid curve) and a case with B = 200 MeV fm -3 
(thin solid curve), in which the maximum gravitational 
mass of a cold neutron star is about 1.85 M Q . For the 
latter case, the metastability timescales lie in a narrow 
range 40-45 s. These, and the metastability masses, are 
both larger than obtained for B = 150 MeV fm~ 3 and 
have narrower ranges. Further increases in B diminish 
the size of the instability window, because the appearence 
of quarks is shifted to progressively larger densities. 




FIG. 4. Lifetimes of metastable stars versus the PNS Mb 
for various assumed compositions. Thick lines denote cases in 
which the maximum masses of cold, catalyzed stars are near 
Ma — 1-45 Mq, which minimizes the metastability lifetimes. 
The thin lines for the npQ and npH cases are for EOSs with 
larger maximum masses (Mg = 1-85 and 1.55 Mq, respec- 
tively. 

Figure |] also compares the metastability time-Ms 
relation found for matter containing hyperons (npH, 
dashed lines fio|| ) or matter with kaons (npK, dotted 
line [pJj) instead of quarks. All three types of strange 
matter are suppressed by trapped neutrinos but 
hyperons always exist in npH matter at finite temper- 
atures and the transition to quark matter can occur at 
lower densities than that for very optimistic kaon cases 
p!l| . Thus, the metastability timescales for npH matter 
can be very short, and those for npK matter are gen- 
erally larger than for npQ matter. Note the relatively 
steep dependence of the metastability time with Mb for 
npH stars, which decreases to very small values near the 
maximum mass limit of hot, lepton-rich, stars. The thick 
npH and npQ lines, as well as the npK line, represent 
minimum metastability times for a given Mb as discussed 
above. The thin npQ and npH lines are for EOSs with 
larger cold, catalyzed maximum mass. 
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Clearly, the observation of a single case of metasta- 
bility, and the determination of the metastability time 
alone, will not necessarily permit one to distinguish 
among the various possibilities. Only if the metastabil- 
ity time is less than 10-15 s, could one decide on this 
basis that the star's composition was that of npH mat- 
ter. However, as in the case of SN 1987A, independent 
estimates of Mg might be available (^0|. In addition, 
the observation of two or more metastable neutron stars 
might permit one to differentiate among these models, 
but given the estimated rate of Galactic supernova (1 
per 30-50 years plj ), this may prove time-consuming. 

Our study has focused on times longer than approxi- 
mately 1 s after core bounce, after which effects of dy- 
namics and accretion become unimportant. Studies of 
the v signal during the first second, during which approx- 
imately 1/3 of the energy is emitted, and at late times, 
as the star becomes optically thin to neutrinos, requires 
more accurate techniques for v— transport. In addition, 
the earliest time periods require the incorporation of hy- 
drodynamics p2[ . 

Our conclusions are that (1) the metastability and sub- 
sequent collapse to a black hole of a PNS containing 
quark matter, or other types of matter including hy- 
perons or a Bose condensate, are observable in current 
and planned v detectors, and (2) discriminating among 
these compositions may require more than one such ob- 
servation. This highlights the need for breakthroughs in 
lattice simulations of QCD at finite baryon density to un- 
ambiguously determine the EOS of high density matter. 
In the meantime, intriguing possible extensions of PNS 
simulations with npQ matter include the consideration of 
heterogenoeus structures [ p3[ , quark matter superfluidity 
p4| and coherent v— scattering on droplets p5| . 
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